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ABSTRACT 


The  generation  of  spontaneous  magnetic  fields  in  a 
plasma  produced  by  a  5.0  Joule,  22  nanosecond  laser  pulse 
and  the  generation  of  the  reversed  magnetic  fields  after 
several  hundred  nanoseconds  due  to  the  interaction  between 
the  laser  produced  plasma  and  a  background  plasma  was 
investigated. 

The  generation  of  the  reversed  magnetic  fields  coincides 
with  the  steepening  of  the  expanding  plasma  front.   This 
steepening  is  caused  by  a  collisional  shock  wave  which  is 
formed  during  the  interaction  of  the  expanding  laser- 
produced  plasma  with  the  photoionized  background  gas.   The 
analysis  of  the  source  term  for  magnetic  field  production 
indicates  that  the  generation  of  the  reverse  fields  is  due 
to  an  axial  temperature  gradient  and  a  radial  density 
gradient.   This  axial  temperature  gradient  is  a  result  of 
the  heating  of  the  plasma  in  the  collisional  shock  wave. 

The  width  of  the  shock  was  measured  to  be  0.2  cm.  The 
magnitude  of  the  maximum  reverse  field  was  measured  at  5.0 
Torr  (H«)  background  pressure  to  be  150  Gauss. 
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I.   INTRODUCTION 

A  dense,  energetic  plasma  can  be  produced  by  focusing 
the  pulse  from  a  powerful  laser  onto  a  solid  target.   Asso- 
ciated with  this  laser-produced  plasma  (hereafter  referred 
to  as  laser  plasma)  spontaneously  generated  magnetic  fields 
have  been  reported  and  investigated  [Refs.  22,24].   These 
investigations  imply  that  spontaneously  generated  magnetic 
fields  arise  generally  during  the  production  of  a  laser 
plasma.   The  strength  of  these  fields  may  vary  with  target 
geometry  [Ref.  31],  background  gas  pressure  [Ref.  3]  and 
amount  of  energy  delivered  to  the  target  by  the  laser.   In 
a  computer  simulation  experiment,  Widner  [Ref.  30]  has 
observed  a  spontaneous  field  intensity  of  the  order  of  1 
Megagauss . 

Stamper  et.  al.  [Ref.  25]  have  reported  spontaneous 
magnetic  fields,  of  the  order  of  0.6  xMegagauss ,  generated 
in  the  plasma  near  the  focus  of  a  high-irradiance  Nd-laser 
beam.   They  used  a  Raman-shifted,  three  channel  Faraday- 
rotation  diagnostic  system  and  observed  that  the  spontan- 
eously produced  magnetic  fields  depend  on  a  variety  of 
experimental  parameters,  such  as  target  material  and 
geometry,  laser  irradiance  and  polarization,  focal  position, 
and  prepulse  effect. 

The  influence  of  magnetic  fields  on  laser  plasmas  is  an 
important  topic  in  laser-fusion  research.   The  presence  of 
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such  fields  in  a  plasma  has  effects  on  electronic  heat 
conduction  and  the  generation  of  plasma  instabilities. 

After  laser  shut-off,  generation  of  reverse  fields  has 
been  observed  at  the  front  of  the  expanding  plasma  for 
ambient  pressures  above  a  certain  critical  value  [Refs.  3 
and  21] . 

This  project  contains  the  study  of  field  reversal  in 
an  expanding  laser  plasma. 

The  spontaneous  and  reversed  fields  were  detected  using 
magnetic  probes  (inductive  loops) .   The  probe  coils  con- 
sisted of  five  turns  of  #4  0  Formvar  copper  wire. 

The  plasma  density  profiles  were  diagnosed  with  elec- 
tric double  probes.   The  double  probes  were  maintained  at 
a  fixed  bias  of  sufficient  magnitude  (20  volts)  to  insure 
collection  of  saturation  ion  current.   Since  the  expansion 
velocity  of  the  laser  plasma  was  greater  than  the  ion 
thermal  velocity,  the  probe  current  was  proportional  to 
the  product  of  the  ion  density  and  the  plasma  expansion 
velocity.   The  expansion  velocity  was  obtainable  from  the 
probe  signals  so  that  the  relative  density  of  the  plasma 
could  be  calculated. 

This  thesis  is  divided  into  7  more  sections.   Section  II 
contains  the  previous  experimental  work  connected  with  the 
field  reversal.   Section  III  discusses  the  basic  theory  of 
spontaneous  magnetic  fields  in  laser  produced  plasmas,  field 
reversal  and  shock  waves  in  a  plasma.   Section  IV  contains 
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the  details  of  the  experimental  methods.   The  experimental 
results  are  presented  in  Section  V.   In  Section  VI,  the 
experimental  results  are  compared  with  the  theory.   Section 
VII  contains  the  summary  and  conclusions.   Recommendations 
for  further  study  are  contained  in  Section  VIII. 


14 


II.   PREVIOUS  RESEARCH 

During  the  latter  part  of  1970,  as  part  of  a  study  of 
the  interaction  between  a  laser  plasma  and  the  ambient 
background  plasma,  Stamper  et.  al.  [Ref.  24]  discovered 
that  magnetic  fields  were  spontaneously  generated  in  the 
laser  plasma.   They  explained  the  generation  of  these 
fields  in  terms  of  thermoelectric  currents  associated  with 
large  temperature  and  density  gradients  near  the  target. 

The  spontaneous  magnetic  fields  of  the  order  of  a  kilo- 
gauss,  were  observed  with  magnetic  probes.   After  this  dis- 
covery, McKee  [Ref.  18]  made  extensive  studies  of  such 
fields  at  the  Naval  Postgraduate  School.   During  McKee ' s 
investigation  the  following  observations  were  made: 

(1)  The  maximum  fields  detected  at  a  fixed  position 
showed  a  systematic  dependence  on  the  ambient  gas  pressure. 

(2)  The  fields  were  symmetric  about  the  target  normal. 

(3)  After  several  hundred  nanoseconds,  a  relatively 
small  magnitude  field  in  a  direction  opposite  to  the  direc- 
tion of  the  initial  field  appeared  behind  the  initial 
fields.   This  reverse  field  appeared  sooner  in  time  for  a 
background  pressure  of  250  m  Torr  than  for  0.1  or  5  m  Torr. 

McKee  developed  a  model  to  explain  spontaneous  magnetic 
fields.   According  to  McKee ' s  model:   If  a  background 
plasma  is  present,  the  expanding  laser  plasma  will  couple 
to  it.   If  the  density  of  the  background  plasma  is  high 
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enough  so  that  the  mean  free  path  for  binary  collisions 
between  the  laser  plasma  ions  and  the  background  plasma 
ions  is  smaller  than  the  dimension  of  the  laser  plasma, 
the  coupling  will  be  collisional  in  nature.   At  lower 
background  pressures  collisionless  coupling  will  arise 
because  the  presence  of  the  spontaneously  generated  magnetic 
field  at  the  front  of  the  laser  plasma  will  lead  to  coupling 
via  the  ion-ion  two  stream  instability  in  the  presence  of 
a  magnetic  field.   As  a  result  of  the  coupling  ambient 
plasma  will  be  swept-up  at  the  front  of  the  expanding 
plasma  and  a  density  shell  will  form.   The  shell  contains 
a  density  gradient  in  the  z-direction  (see  Figure  1) ;  then 
field  generation  can  occur  in  the  density  shell  if  a  radial 
temperature  gradient  arises  (see  Section  III  for  details). 

After  McKee,  Bird  developed  another  model  (interaction 
model)  to  explain  the  pressure  dependence  of  spontaneous 
magnetic  fields.   This  model  attributes  the  increase  in  the 
magnitude  of  the  initial  spontaneous  magnetic  fields  with 
increasing  background  gas  pressure,  to  the  interaction  which 
occurs  between  the  laser  plasma  and  the  ambient  background 
plasma.   The  interaction  increases  the  axial  density  gradi- 
ents of  the  front  of  the  expanding  plasma  and  this,  in  turn, 
increases  the  rate  at  which  spontaneous  magnetic  fields 
are  generated  in  the  interaction  region. 

During  Bird's  investigation  no  density  shell  was  detected 
Also  the  shell  model  could  not  explain  the  field  reversal. 
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The  following  observations  on  reverse  field  generation 
have  been  reported  by  Bird: 

(1)  Field  reversal  is  coincident  with  a  strong  inter- 
action between  the  laser  plasma  and  the  background  plasma 
as  is  evidenced  by  a  steepening  and  slowing  of  the  front 
of  the  expanding  plasma. 

(2)  The  onset  of  field  reversal  occurs  at  the  front 
of  the  expanding  laser  plasma  and  it  occurs  earlier  in 
space  and  time  at  higher  pressures. 

(3)  Reverse  field  can  be  produced  long  after  laser 
shut-off  if  the  laser  plasma  encounters  a  barrier. 

(4)  Reverse  field  is  initially  produced  at  the  front 
of  the  laser  plasma  streaming  into  a  background  of  H2  at 
pressures  above  450  m  Torr  before  any  significant  pile-up 
of  the  laser  plasma  has  occurred.   This  narrow  component 
of  reverse  field  appeared  to  be  generated  in  the  swept-up 
hydrogen  background  plasma  and  it  was  followed  by  the 
appearance  of  a  broad  component  of  reverse  field  as  the 
plasma  front  steepens  and  slows  down  due  to  the  interaction 
between  the  laser  plasma  piston  and  the  background. 

These  results  suggest  that  the  generation  of  the  reverse 
magnetic  fields  can  be  explained  by  formation  of  a  shock 
wave  driven  by  the  expanding  laser  plasma  piston.   Such  a 
shock  wave  compresses  and  heats  the  plasma.   The  formation 
and  the  effects  of  the  shock  waves  in  a  plasma  are  discussed 
in  the  following  section. 
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III.   THEORY 

A.   LASER  PRODUCED  PLASMAS 

The  production  of  plasma  by  irradiation  of  a  solid 
target  was  first  treated  theoretically  by  Basov  and  Krokhin 
[Ref.  1] . 

The  plasma  is  created  by  heating,  vaporization  and 
ionization  of  the  target  material.   The  overdense  plasma 
(n   >  n  )   at  the  surface  of  the  target  expands  until  the 
electron  density  in  the  front  drops  below  the  critical 
density  n  •   Strong  absorption  occurs  when  the  density  is 
near  the  critical  density  n   and  the  plasma  is  heated  to  a 
high  temperature  and  the  pressure  gradient  accelerates  the 
plasma.   After  the  laser  shut-off,  the  plasma  front  expands 
at  a  constant  velocity.   Thus,  three  distinct  states  can  be 
observed  during  the  course  of  the  expansion  of  the  laser 
plasma. 

When  the  expansion  is  into  an  ambient  background  gas 
rather  than  into  vacuum,  the  last  stage  looks  quite  different 


The  critical  density  n   for  the  laser  plasma  is  the 
density  at  which  the  electron  plasma  frequency  is  equal  to 
the  frequency  of  laser  radiation  (n   =  oo  mge  /e^)  .   At  elec- 
tron densities  above  this  critical  value,  the  laser  plasma 
reflects  most  of  the  incident  radiation.   The  critical 
density  for  1.06  ym  radiation  is  about  10    cm~^. 


At  this  stage,  for  a  background  gas  pressure  greater  than 
0.2  Torr,  Bobin  et.  al.  [Ref.  2]  observed  a  luminous  sheet 
ahead  of  the  expanding  plasma.   This  luminous  sheet  lasted 
more  than  150  nsec  after  the  laser  pulse  peak  value.   Its 
velocity  decreased  with  time.   The  authors  attributed  the 
origin  of  the  luminous  sheet  to  a  shock  wave  driven  by  the 
plasma  acting  as  a  piston,  because  their  observed  time  evolu- 
tion of  the  luminous  front  agreed  both  qualitatively  and 
quantitatively  with  a  simple  shock  wave  model. 

A  one-dimensional,  simplified  model  treating  the  expan- 
sion of  the  laser  plasma  into  vacuum  has  been  given  by 
Ready  [Ref.  20].   In  this  treatment  it  is  assumed  that: 

(1)  The  plasma  expands  in  the  z-direction. 

(2)  Incident  radiation  is  absorbed  via  the  inverse 
Bremsstrahlung  process  (free-free  absorption  by  electrons 
in  the  laser  plasma. 

(3)  The  target  is  initially  a  cold  fully  ionized 
plasma. 

(4)  Radiation  losses  can  be  neglected. 

(5)  The  electrons  and  ions  have  a  common  temperature 

(6)  Temperature  is  independent  of  position. 


2 

The  equipartition  time  for  transfer  of  absorbed  energy 

from  electrons  to  ions  can  be  approximated  by 

Tf-i  *25(-i)Te3/2(°K)  [n  (cnT^Z]'1 
ei    m  c 

where  m.  is  the  ion  mass  and  m   is  the  proton  mass.   For 

4-  P 

n  =  1021  and  Te  =  100  ev,  t£.  <  1  nsec. 
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(7)   The  plasma  expansion  velocity  is  of  the  form 


dtV 


pysically,  the  pressure  gradients  within  the  laser  plasma 
are  driving  the  expansion.   The  form  of  the  pressure 
variation  can  be  approximated  by 


P   =   Ps(l  -  §)        Z  <  Z 


where  P   is  the  pressure  near  the  surface  of  the  target, 

and  Z  is  the  distance  of  the  plasma  front  from  the  target 

surface . 

The  temperature  of  the  plasma  can  be  determined  by 
i 

equating  the  rate  at  which  the  laser  beam  deposits  energy 
in  the  plasma,  to  the  rate  of  change  of  plasma  thermal 
energy  plus  the  rate  at  which  work  W  is  done  by  the  plasma 
For  one  dinensional  plasma  of  area  A,  expanding  in  the  z- 
direction  this  gives : 


1  dW      Id,      N1  m         3P    dz  ,,  . 

A  dz      2  dt   e  'i'      3z    dt 


From  the  plasma  equation  of  state 


P   =   (ni  +  ne)kT  (2 


one  can  obtain  the  variations  of  electron  and  ion  densities 
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ni   =   nis(1  -  f}  ;   ne   =   nes(1  "  l}        (3 


where  n  •   and  n.    are  the  ion  and  electron  densities  at  the 

is      es 

target  surface.   The  work  done  by  the  plasma  in  the  expan- 
sion increases  the  plasma  kinetic  energy  so  that: 


8P   dz      d  ,1,     ,     ,  /dz,  2,  ... 

li"  '  dt   =   dt[2(riimi+rieme)  (dt)  ]  (4) 


By  integrating  Equations  (1)  and  (4)  from  zero  to  Z  and 

the  result, 

Z 
N  =  A  /  ns(l  -  |)dz  =  Ansz/2  , 
0 

where  N  is  the  total  number  of  ions  or  electrons  in  the 

plasma  and  n   is  the  corresponding  number  density  at  the 

target  surface,  Ready  obtained  the  equations 


(N.m. +N  m  ) Z     ,2^ 

T   =   [   x  x   e  e   ]  •  14  (5) 

L  6K(Ni+Ne)   J    dt2  lDj 


,3  2 

2-^-  =   [24(Nimi+Neme) ]W  (6) 


These  two  coupled  equations  were  solved  numerically  by 
McKee  for  the  conditions  of  this  experiment.   The  maximum 
temperature  obtained  was  of  the  order  of  100  eV  and  this 
maximum  occurred  a  few  nanoseconds  after  the  arrival  of  the 
maximum  laser  intensity  at  the  target. 
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In  the  above  treatment,  it  was  assumed  that  the 
expansion  of  the  plasma  is  into  a  vacuum.   The  presence 
of  a  background  gas  (pressures  above  1  m  Torr)  influences 
the  expansion  dynamics. 

B.   INTERACTION  BETWEEN  COUNTERS TREAMING  PLASMAS 

The  optical  radiation  from  a  dense,  hot  laser  plasma 
is  known  to  be  rich  in  the  far  ultraviolet  and  near  x-ray 
regions  of  spectrum.   Line  radiation  and  Bremsstrahlung 
radiation  from  a  plasma  of  T  =  100  eV  is  energetically 
capable  of  ionizing  an  ambient  background  gas. 

Koopman  [Ref.  16]  has  investigated  the  photoionization 
produced  from  a  plane  copper  target  by  an  8-Joule,  20  nsec 
laser  pulse  at  6943  A°  from  a  Ruby  laser.   He  found  that 
the  ambient  gas  was  almost  completely  ionized  to  distances 
of  the  order  of  a  centimeter  for  a  background  gas  pressure 

of  the  order  of  0.2  Torr  and  at  large  radii  the  degree  of 

-2 
photoionization  fell  as  r   .   Total  optical  energy  emitted 

between  1000-2000  A°  was  estimated  as  approximately  1±0.5 
Joule.   Therefore,  when  a  laser  plasma  expands  into  a 
photoionized  ambient  background  gas,  plasma-plasma  inter- 
action must  be  considered  on  the  expansion  dynamics. 

According  to  Dean  [Ref.  7]  during  the  period  of  the 
laser  irradiation  of  the  target  (early  time) ,  the  dynamics 
of  the  laser  plasma  expansion  corresponds  to  the  propagation 
of  a  radiation-driven  detonation  wave.   After  the  end  of 
the  laser  pulse  (late  time)  the  plasma  expands  as  a  blast 
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wave.   This  model  assumes  that  at  early  times  of  the 
expansion  phase,  there  is  momentum  transfer  between  the 
expanding  front  and  the  ambient  background  plasma  and  that 
as  the  front  propagates  into  the  background  plasma,  further 
energy  from  the  laser  is  absorbed  behind  the  front  in  which 
the  density  is  near  the  critical  density  and  is  transported 
to  the  front  by  electron  heat  conduction.   After  the  end 
of  the  laser  pulse,  all  the  momentum  in  the  driving  laser 
plasma  piston  is  assumed  transferred  to  the  swept-up  back- 
ground plasma  which  propagates  as  a  blast  wave. 

Hall  [Ref.  12]  investigated  the  dynamics  of  the  expan- 
sion of  the  laser  plasma  into  an  ambient  background  of 
argon.   During  his  investigation,  a  0.5  Joule,  20  nsec  Ruby 
laser  was  used.   The  target  was  a  Tantalum  slab.   Using  an 
image-converter  camera  and  double  probes  for  diagnostics, 
he  found  that  the  late  time  expansion  dynamics  of  the  laser 
plasma  could  be  explained  by  blast  wave  theory. 

Koopman  [Ref.  16]  has  suggested  that  the  observed  momen- 
tum coupling  is  the  result  of  multi-coulomb  collisions 
occurring  between  particles  of  the  two  plasmas.   If  a  strong 
momentum  transfer  interaction  occurs,  one  would  expect 
"snowplowing"  or  displacement  of  the  ambient  plasma  by  the 
leading  edge  of  the  laser  plasma  piston,  with  a  corres- 
ponding pile-up  of  piston  plasma  behind  the  interface  between 
the  laser  plasma  and  the  background  plasma.   The  ambient 
and  piston  plasmas  would  be  intermixed  in  the  region  of 
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interface,  that  is,  the  model  is  one  in  which  the  background 
plasma  is  "captured"  at  the  front  of  the  expanding  plasma, 
increasing  the  total  mass  of  the  streaming  plasma  and,  by 
conservation  of  momentum  decreasing  the  expansion  velocity. 

A  one  dimensional  expansion  model  can  be  justified  on 
the  basis  of  the  following  considerations.   At  the  target 
surface  the  laser  beam  diameter  is  about  1.5  mm  and  the 
laser  radiation  is  absorbed  within  a  distance  of  the  order 
of  0.01-0.1  mm  from  the  target  [Ref.  20].   Therefore,  the 
ablated  target  material  initially  will  have  larger  dimen- 
sions in  the  radial  direction  than  in  the  axial  direction. 
Since  the  expansion  of  the  laser  plasma  is  driven  by  pressure 
gradients; 


(VP)  |   -   f     and     |(VP)  |   -   f 

z  r 


from  the  above  discussion,  5   >>  5  ,  therefore 

r     z 

I  (VP)    >>  I  (VP)  I  where  5   and  6   are  characteristic 

I       z  I        I       r  I  r  z 

radial  and  axial  scale  distances.   From  the  above  inequality, 
a  one  dimensional  laser  plasma  expansion  can  be  approximated. 
When  the  laser  plasma  expands  into  an  ambient  background 
plasma,  the  axial  position  of  the  expanding  front  will  be 
governed  by  the  interaction  between  the  two  plasmas.   An 
empirical  expression  for  the  late  time  (after  laser  shut- 
off)  dependence  of  the  axial  position  of  the  plasma  front 
on  time  and  on  the  ambient  (nitrogen)  plasma  density,  p  , 
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is  given  by  Bird  [Ref.  3]: 


_  ,  1  . 0.10.5  ,-» 

Za(— )    t  (7) 


C.   INITIAL  PRODUCTION  OF  SPONTANEOUS  MAGNETIC  FIELDS 

In  the  presence  of  electric  and  magnetic  fields,  the 
equation  of  motion  for  the  electron  fluid  is  given  by 
generalized  Ohm's  law. 


8V 

n  m  [ttt~+  (V  -V)V  ]   =   -en  [E  +  V  XBJ-V-P  +  P  .    ( 
e  e  9t     e    e         e     e       e    ei 


where  V  (r,t)  is  the  fluid  velocity  of  electrons,  E(r,t) 

and  B(r,t)  are  the  averaged  electric  and  magnetic  fields, 

P   is  the  kinetic  stress  tensor  and  P    represents  the 
e  ei 

colLisional  momentum  transfer  between  electrons  and  ions 

per  unit  volume  and  unit  time. 

If  one  assumes  that  the  electron-electron  collision  time 

t    is  much  less  than  the  characteristic  time  for  field 
ee 

production  which  is  of  the  order  of  duration  of  the  laser 

pulse,  then  V*P   reduces  to  the  electron  pressure  gradient, 

VP  .   If  the  macroscopic  flow  velocity  of  the  electron 

fluid  element  is  small  compared  to  the  random  thermal 

velocity  of  the  electrons  then  the  (V  -7)V   term  can  be 
1  e    e 

neglected.   For  the  aluminum  laser  plasma,  the  expansion 

7 
velocity  of  the  front  is  of  the  order  of  -10   cm/sec 

[Ref.  5] .   If  one  assumes  an  electron  temperature  of  100  eV, 
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then 


|V  | 

t-jj      ~   0.01  . 

[kT  /m  ]1/z 
e   e 


Therefore  the  above  approximation  is  justified. 

If  the  frequencies  of  interest  are  much  less  than  the 
electron  plasma  frequency,  i.e.,  co  <<  go    the  electrons 
can  be  considered  massless  and  react  instantaneously  to 
applied  fields. 

P  .  can  be  approximated  by  the  resistive  drag  term, 

n  eJ  el 

e  .  . 

— - — ,  where  a  is  the  conductivity  of  plasma.   Finally, 

eq.  (8)  can  be  written 


E   =   -  -  V  x  B — VP  (9) 

a    e      nee 
e 


By  taking  the  curl  of  E  and  using  the  Maxwell  equations, 
eq.  (9)  becomes 


-VxE   =   4^   =   -  -Vxj  +  V  X(V  xB)  +  V  x  — VP  (10 

dt       a         e       -  n  e   e 

e 


Also,  from  the  Maxwell  equations, 


VXB   =   UQJ  +  -^  $f  (ID 

c 


The  ratio  of  the  two  terms  on  the  right  hand  side  of 
eq.  (11)  is  of  the  order 
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1 

• 

2 

E 

c 

^o 

J 

1_  _E_ 

2     T 


I?)  (12) 


-  \x    oE  '  v     3t 

o 

_p 

where  t   is  the  duration  of  the  laser  pulse  (t  *  5  x  10    sec) 
o  o 

for  all  situations  of  interest  a    >  1  mho.   Therefore, 
eq.  (12)  becomes: 


<   f°   <   X0^   <     -4  <K      ± 


2  QT         T  T 

c  y  ax        o      o       o 
o  o 


The  above  inequality  shows  that 


_L   II 
2  '  at 


^oJ 


i.e.,  the  displacement  current  can  be  neglected  in  compari- 
son to  the  conduction  current.   In  doing  so,  equation  (11) 
becomes 


V  x  b   =   u   J  (13) 


using  eq.  (13),  the  first  term  in  eq.  (10)  becomes 


-  i  Vxj   = L_  Vx(7xB)   = — [V(7«B  -  v2B)  ] 

a  u  a  u  a 

o  o 

(13) 


since 
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7-  B   =   0  ,     -  -|vx  J   =   -i~  V2  B        (14) 


The  third  term  in  eq.  (10)  can  be  written 


Vx-^VP         =    V    x   ~=-V(n    kT    )       =    V    x(-7T     + — — Vn    ) 
nee  n   e        e      e  eenee 

e  e  e 

kT 
=  v  x(F-|vne)  (15) 

vne 

The  density  gradient  is  a  vector;  =  A  but  the  tempera- 

kT     ne 
ture  is  a  scalar  quantity;  =  0   then,  using  the  vector 


identity; 


Vx(0A)   =   j2(YxA  +  (V0)xA  (16) 


Eq.  (15)  becomes: 
kT 


Vx( — ^Vn  )   =   -VT  x-i-Vn  (17) 

nee      e   e  n    e 
e  e 


putting  all  terms  in  their  final  form  together,  eq.  (10) 
becomes 


If   =   tAv2b  +  Vx(veXB)  +  ^!VTe)x(;r-vV]     (18) 

o  e 


The  first  two  terms  on  the  right  hand  side  of  eq.  (18) 
are  the  diffusion  and  convection  terms.   There  is  no  external 
magnetic  field  applied,  i.e.,  initially  B  =  0.   Therefore 
the  generation  of  an  initial  magnetic  field  requires  that 
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the  last  term,  the  thermal  source  term  (S)  be  non-zero 
due  to  non-parallel  temperature  and  density  gradients. 

The  initial  configuration  of  magnetic  field  will  decay 
away  in  a  diffusion  time  t  [Ref.  14], 


x  =  iisti  d9) 


where  L  is  a  length,  characteristic  of  the  spatial  varia- 
tion of  B.   For  times  short  compared  to  the  diffusion  time 
or,  in  other  words,  when  the  conductivity  is  so  large  that 
the  first  term  in  eq.  (18)  can  be  neglecged,  equation  (18) 
can  be  reduced  to: 


s  VkT      Vn 

|f  =  V  x(VexB)  +  (— g-£)  x  (_®)  (20) 


This  equation  expresses  the  fact  that  the  rate  at  which  the 
fields  vary  in  a  region,  depends  on  the  rate  at  which  the 
fields  are  created  in  the  region  and  the  rate  at  which  the 
field  flows  in  or  out  of  the  region. 

For  adiabatic  expansion  of  the  electrons,  PV   =  const. 

M 
e   y 

or  n  kT  ( )   =  const,  (where  M   is  the  total  electron 

e   e  m  n  e 

e  e  , 

1— Y 

mass  in  a  given  volume  V)  or  kT  n      =  const,  and  the 

e  e 

source  term  vanishes  (VT  *vn  aVn  *7n   =  0) .   Thus,  non- 

e    e    e    e 

adiabatic  conditions  are  necessary  for  the  generation  of 
magnetic  fields  by  the  thermal  source  term. 

If  the  laser  energy  is  supplied  to  the  electrons  at  a 
sufficiently  high  rate,  then  the  magnetic  source  term  is 
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non-zero  in  the  resulting  non-adiabatic  state.   Stamper 
[Ref.  26]  analyzed  the  problem  by  means  of  the  electron 
heat  equation. 


dWH   -    *-  •  S  -  k  T  33,  (21) 


dt      y-1    dt         dt 


where  Wu  is  the  energy  supplied  to  the  electrons  per  unit 

n 

time.   The  first  term  on  the  right  hand  side  of  equation  (21' 

nk    dT   dWI 

— —t-    '    j—  =  -rrr  where   WT  is  the  internal  energy  and  the 

y-1    dt    dt  I  3J 

second  term  on  the  right  hand  side  of  equation  (21) 

d     dWW 
kTS-z-  =  -T-—     where  WTT  is  the  work  done  by  the  electrons, 
dt    dt  W 

According  to  Stamper,  the  ratio  W  /W   is  the  measure  of 

n   I 

non-adiabatic  conditions.   Stamper  approximated  Wg  by  th 
product  of  the  absorption  coefficient  k  (cm-1)  and  the 

radiation  intensity  I  (Watt/cm  ) ,  while  Wj  was  approximated 

nekT 

by  __ where  t   is  an  average  electron  heating  time.   If 

xh  h 

one  assumes  kv  is  due  to  inverse  Bremsstrahlung  absorption, 

in  the  underdense  region,  that  gives: 

W      t  e.     ,  — 
H      ei    I/cn  ,99» 

"  t  .  '  n  kT  [ZZ) 

W       ei    e 

4  —  . 
where  T    is  the  electron-ion  collision  time,  r\    is  the 
ei 


4        2.7  x  io  2    _  3/2  ,0_.  -         ,n21    -3 

t    =  5 •  Te    (°K)  for  n  ■  ~  10    cm 

ei      ,   -3.  i 

ni(cm   ) 

-14 
and  T   -  100  eV.   From  footnote  (2),  x  .  -    3  x  10     sec. 

ID  ei  2 

t£  .  -  8.1  x  10     sec,  and  the  ratio  t£./t  .  ~    9.25  x  10   -—  >>  1 
ei  ei   ei  -  m 
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index  of  refraction  of  plasma  (r\    ~  1)  .   Actually  near  the 
cut-off,  the  index  of  refraction  goes  to  zero,  thus  the 
ratio  W^/W.  >>  1  and  the  plasma  appears  non-adiabatic . 
During  the  time  the  plasma  is  absorbing  the  laser  radiation, 
for  the  conditions  of  this  investigation,  the  ratio  W^/W 

rl   J. 

is  of  the  order  of  one  and  the  plasma  appears  non-adiabatic. 

•  •  • 

But  if  W   <<  W  and  W„  then  the  plasma  appears  adiabatic. 
The  direction  of  magnetic  field  is  determined  by  the 
geometry  of  the  laser  plasma.   At  early  times,  the  laser 
plasma,  as  justified  before,  will  expand  in  the  axial  direc- 
tion and  create  an  axial  density  gradient.   The  temperature 
gradient  is  due  to  the  finite  radial  extent  of  the  laser 
beam  and  is  in  the  negative  radial  direction.   The  cross 
product  of  these  two  gradients  will  generate  azimuthal  and 
axisymmetric  magnetic  fields  [Fig.  6]  as  observed. 

D.   FIELD  REVERSAL 

At  later  times,  after  cessation  of  the  laser  pulse, 
the  gradients  become  smaller  and  field  production  becomes 
insignificant.   However  if  the  laser  plasma  expands  into 
a  background  plasma  of  sufficient  density,  shock  heating 
at  the  front  will  produce  an  axial  temperature  gradient 
[Ref .  21] .   Since  a  radial  density  gradient  always  exists 
in  the  plume-shaped  laser  plasma,  magnetic  field  in  a 
direction  opposite  to  the  initial  field  will  be  generated. 
The  equation  describing  this  field  is: 


kT      V  n 
£|  =   (V   _£)  x  {JL^  (23 

<3t        z   e         n 
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Schwirzke  [Ref.  21]  and  Bird  [Ref.  3]  observed  this  magnetic 
field  reversal  several  hundred  nanoseconds  after  laser  shut- 
off.   The  production  of  reverse  field  at  the  front  of  the 
laser  plasma  requires  that  the  ratio: 

I  v  t  I  I  v  n 

1  z  e ' '  re' 


V  T    V  n 
re11  z  e 


be  greater  than  one. 

The  gradients  in  this  ratio  can  be  approximated  by: 


T  ri  T  n 

V  T   =  ,3-  ,    V  n   =  tS-  ,    7  T  z    n3-  ,    V  n 


z  e    6  _  '     z'e  ~  6    '     re  "    6  _  '     r  "e  ~  5  n 
zT  zn  rT  rn, 


where  the  5's  are  the  characteristic  gradient  scale  lengths. 

The  ratio  becomes  6   *5  _/6  _*5    .   For  a  one  dimensional 

zt]      rT   zT   m 

adiabatic  compression  or  shock  compression  in  the  z-direction, 

5-^6    and  5  _  <<  5  _.   The  large  electron  heat  conduc- 
zT     zt]  zT      rT  3 

tion  provides  that  6  _  >  6    and  thus  the  ratio  becomes 
greater  than  one  ie. ,  reverse  field  production  is  possible 
under  these  conditions. 

The  production  of  reverse  field  appears  to  be  a  conse- 
quence of  the  interaction  between  the  laser  plasma  and  the 
ambient  background  plasma.   The  interaction  results  in 
slowing  down  and  a  steepening  of  the  front  of  the  laser 
plasma  (pile-up) . 

Equation  (23)  predicts  that  the  field  reversal  will  occur 
sooner  and  grow  to  a  larger  magnitude  as  the  background 
pressure  increases.   Bird  [Ref.  5]  justified  this  prediction 
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and  also  observed  that 

(1)  The  mechanism  for  production  of  reverse  fields 
is  "turned-on"  sooner  at  higher  background  pressures. 

(2)  Field  reversal  is  coincident  with  a  strong  inter- 
action between  the  laser  plasma  and  the  ambient  background 
as  is  evidenced  by  a  steepening  and  slowing  of  the  front 
of  the  expanding  plasma. 

The  electron  temperature  required  to  produce  a  given 
field  reversal  can  easily  be  estimated  by  dimensional 
analysis  of  eq.  (23) .   From  that  equation 


K  T       BQ     *  5   •  6 

e   ^    8max r z  (?&\ 


where  6   and  6   are  radial  and  axial  characteristic  scale 
r      z 

lengths  and  t  is  the  characteristic  time  for  reverse 
magnetic  filed  production.   Bird  estimated  the  required 
temperature  to  be  of  the  order  of  1  eV  for  the  parameters 
of  his  investigation  which  were  similar  to  the  ones  of  the 
present  investigation. 

E.   SHOCK  WAVES  IN  A  PLASMA 

If  a  large  concentration  of  energy  in  a  continuous 
medium  is  suddenly  released,  it  will  spread  into  the  medium 
and  at  its  forefront  will  be  a  shock  wave.   A  shock  wave  is 
a  sudden  but  continuous  change  in  the  state  of  medium,  heats, 
compresses  and  accelerates  the  medium  through  which  it 
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propagates.   The  shock  wave  may  propagate  through  a 
plasma. 

The  normal  process  of  shock  formation  is  the  pressure 
dependence  of  wave  velocity,  so  that  a  finite  amplitude 
disturbance  with  increasing  pressure  steepens  as  the  higher 
pressure  portions  move  faster  than  the  smaller  ones  in 
front.   In  order  to  obtain  a  steady  profile  wave,  some 
mechanism  must  exist  to  limit  the  steepening.   This  limi- 
tation on  steepening  can  be  provided  by  either  dissipation 
or  dispersion.   A  steady  state  shock  occurs  when  a  diffu- 
sion process  balances  the  steepening.   Jump  conditions  are 
then  obtainable  from  the  conservation  of  mass,  momentum 
and  energy.   In  gas  dynamics,  the  flows  are  properly  des- 
cribed by  the  Navier-Stokes  equations.   Any  shock  wave 
present  in  the  flow  must  therefore  also  be  contained  in  the 
solutions  of  the  Navier-Stokes  equations,  under  appropriate 
initial  and  boundary  conditions  [Ref.  3].   However  if  the 
thickness  of  the  shock  wave  becomes  so  small  as  to  be 
comparable  to  the  mean  free  path,  the  Navier-Stokes  equations 
become  invalid  in  describing  the  flow  in  the  shock.   One 
must  then  use  the  Boltzmann  equation  or  similar  kinetic 
equations  for  the  description  of  the  shock  structure. 

The  basic  features  of  the  shock  structure  in  a  plasma 
are  related  to  the  slow  character  of  the  energy  exchange 
between  ions  and  electrons  and  to  the  high  electron  mobility, 
as  a  result  of  which  the  electron  heat  conduction  greatly 


34 


exceeds  the  ion  heat  conduction.   Maxwellian  distributions 
in  the  electron  and  ion  fluids  are  established  quite 
rapidly,  in  a  time  of  the  order  of  the  time  between  parti- 
cle collisions  [Ref.  33].   On  the  other  hand  the  equili- 
bration of  the  temperatures  of  ion  and  electron  fluids 
takes  place  much  more  slowly,  because  of  the  large  differ- 
ence between  the  electron  and  ion  masses.   This  relaxation 
process  determines  the  shock  front  thickness  in  a  plasma. 

A  qualitative  discussion  of  the  consequences  of  the 
low  rate  of  energy  exchange  between  electrons  and  ions  is 
given  by  Zel'dovich  and  Raizer  [Ref.  33].   In  this  simpli- 
fied discussion  it  is  assumed  that: 

(1)  electron  heat  conduction  does  not  differ  from  ion 
heat  conduction 

(2)  ionization  does  not  take  place  in  the  shock  .wave; 
the  wave  is  propagated  through  a  gas  which  is  already 
ionized 

In  a  reference  frame  moving  with  the  wave  a  considerable 
part  of  the  kinetic  energy  of  the  gas  entering  the  com- 
pression shock,  is  irreversibly  converted  into  heat  through 
the  action  of  ion  viscous  forces.   The  increase  in  ion 

temperature  across  the  compression  shock  is  of  the  order 

m.V^ 
AT.  ~  — r —  where  V  is  the  velocity  of  the  incident  gas, 

equal  to  the  shock  front  velocity  and  m.  is  the  mass  of  an 

ion.   The  thickness  of  the  viscous  shock  is  determined  by 

the  time  between  ion  collisions,  x . ;  it  is  of  the  order  of 
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an  ion-ion  collision  mean  free  path  & . ~v.,t .  where  v^, 

r     1   th  1        th 

is  the  thermal  velocity  of  the  ions  in  the  compression 

shock.   During  the  compression  time  t.  the  ions  do  not 

succeed  in  transfering  any  appreciable  thermal  energy  to 

the  electrons,  since  the  characteristic  exchange  time 

m.  i  /~ 
t  .  -    ( — )    t .  is  very  large  compared  to  the  ion  collision 
e 

time  x . . 

l 

The  increase  in  electron  temperature  across  the  shock 
due  to  the  conversion  of  the  kinetic  energy  of  the  elec- 
trons into  heat  through  the  action  of  electron  viscous 

forces  is  negligibly  small.   It  is  of  the  order  of 

2 

AT   ~  meV 

e   — : —  ,  that  is  smaller  by  a  factor  of  m  /m.  than  AT. . 
k  2  e   i        l 

The  electron  gas  in  the  compression  shock  is  heated  by  some 
other  process.   Since  the  dissipation  of  energy  due  to 
electron  viscous  force  is  very  small,  Zel'dovich  and  Raizer 
[Ref.  33]  suggested  that  the  compression  and  heating  of 
the  electron  gas  in  the  compression  shock  is  adiabatic: 


T  0       p~    , 

7^   =   i^)1'1  (25 

Tel      pl 


where  T  .  and  T  ~  are  the  electron  temperatures  in  State  I 
el      e2  r 

and  State  II,  p,  and  p2  are  the  electron  densities  in  State 

I  and  State  II  and  y  is  the  ratio  of  specific  heats.   For 

P2 
a  strong  shock  wave,  the  ratio  —  is  of  the  order  of  4.   Thus 

T  pl 

e2 
= —  ~  2.5  for  y  =  1.4.   Therefore  in  a  strong  shock  wave 

el 

propagating  through  a  plasma  with  equal  ion  and  electron 
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temperatures,  behind  the  shock  there  is  considerable 
difference  in  the  temperatures  of  the  two  gases.   After  a 
small  portion  has  undergone  the  shock  compression  it  begins 
to  transfer  thermal  energy  from  the  ions  to  the  electrons 
and  this  leads  to  an  equilibration  of  the  temperatures  in 
a  time  of  the  order  of  the  energy  exchange  time  t  . .   The 
final  temperature  T,  (=  T   =  T.)  is  determined  by  the  general 
equations  of  conservation  for  shock  fronts.   The  thickness 
of  the  relaxation  region  behind  the  shock  in  which  plasma 
approaches  equilibrium  (T   =  T.  =  T, )  is  of  the  order  of 

£  Pl 

AX  -  vxTei       where       v±      =      — V  . 

Jaffrin  and  Probstein  [Ref.  15]  investigated  the  structure 
of  a  one  dimensional  shock  wave  in  a  fully  ionized  plasma 
consisting  of  electrons  and  singly  ionized  ions  in  the 
absence  of  any  applied  external  electric  or  magnetic  fields. 

The  analysis  of  the  shock  structure  is  based  on  a  two- 
fluid  continuum  approach  in  which  they  took  as  a  "model" 
the  Navier-Stokes  equation,  including  the  electric  field  and 
transfer  terms.   In  so  doing  the  overall  width  of  the  shock 
transition  region  turns  out  to  be  considerably  larger  than 
the  characteristic  mean  free  path.   They  found  that 

(1)   for  mach  number  M  <  1.12  (weak  shock)  the  shock 
profile  is  monotonic,  with  the  shock  thickness  of  the  order 
of  ly/e   where  &_  is  the  down  stream  (behind  shock)  mean  free 
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m   .  ,„ 
e  1/2 
path  of  ion-ion  collisions  and  e  =  (— )    .   The  important 

i 

dissipative  mechanisms  are  electron  thermal  conduction  and 

the  energy  exchange  between  ions  and  electrons. 

(2)  M  >  1.12  (moderate  to  strong  shocks).   There  is 
still  a  thermal  shock  layer  whose  approximate  thickness 

is  l~/z.      In  this  layer  the  electron  temperature  increases 
and  approaches  its  downstream  value  while  the  ion  tempera- 
ture and  velocity  changes  are  relatively  small.   However, 
imbedded  in  this  first  layer  there  exists  a  second  inner 
shock  whose  thickness  is  of  the  order  of  the  downstream 
ion-ion  collision  mean  free  path  £„.   In  this  second  layer, 
the  dissipations  are  due  to  the  ion  viscosity  and  thermal 
conduction.   The  heating  and  compression  of  ions  occur 
mainly  in  this  second  shock  layer.   However,  because  of 
the  ions'  much  heavier  mass,  the  ions  do  not  have  time  to 
exchange  energy  with  the  electrons  in  a  few  mean  free  paths 
and  the  electron  temperature  is  essentially  unaffected. 
The  ion  and  electron  temperatures  then  equilibrate  in  the 
downstream  continuation  of  the  thermal  layer.   The  thick- 
ness of  this  downstream  layer  is  of  the  order  of  'l~/e. 

(3)  Debye  length  of  the  order  of  the  mean  free  path 
(A   -  ^2^*   In  tn:>-s  case  spatial  oscillations  of  wave- 
lengths which  are  of  the  order  of  the  Debye  length  are 
present.   The  charge  separation  is  no  longer  small  and  the 
electrical  effects  interact  with  the  flow  in  the  two 
electric  shock  layers . 
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The  structure   of   a   collisional   weak   shock   in 

the   presence  of  magnetic   field  which  is  in  the  plane  of 

shock  was  studied  by  Hu  [Ref.  13].    He   found   that 

the  shock  thickness  depends  on  four  parameters,  the 

electron  mean  free  path  X    ,  the  electron  1 armor  radius  r 

e  e 

based  on  thermal  speed,  the  electron  penetration  length 

m  c^  _ 

e  2 

r  =  * — t— TR-,  and  the  electron-ion  mass  ratio  z    .   Hu 

p    , .     2. 1/2 

e         (47Tri  e  ) 

e 

discovered  two  important  features  of  the  shock  structure. 

(1)  It  is  possible  to  have  collisional  shocks  with 
oscillatory  structure  due  to  electron  inertia. 

(2)  It  is  possible  to  have  collisionless  shock  with 
collisions  whose  thickness  is  smaller  than  the  collision 
mean  free  path  due  to  the  complicated  interactions  between 
electromagnetic  fields  and  the  plasma.   Often  in  interpreting 
experiments,  shocks  exhibiting  these  two  features  are  con- 
sidered as  collision-free  shocks. 

In  general,  the  mean  free  path  for  collisions  in  a 
plasma  determines  the  type  of  the  shock.   A  shock  layer, 
which  depends  on  the  interaction  length,  can  be  thinner 
than  a  mean  free  path.   These  shocks  provide  examples  for 
studying  the  mechanism  of  collisionless  dissipation.   A 
similar  definition  of  the  collisionless  property  was  given 
by  Friedman  [Ref.  10].   For  a  given  shock  structure,  the 
dissipation  provided  by  classical  collisions  (ordinary 
momentum  transfer  collisions)  must  be  less  than  the  dissi- 
pation which  is  required  to  form  a  steady  shock. 

The  first  conclusive  evidence  of  the  collisionless 
shocks  come  from  astrophysics.   The  plasma  flow  from  sun, 
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the  solar  wind,  is  highly  supersonic  when  encountering  the 
magnetic  field  of  the  earth,  and  a  shock  wave  must  be 
formed.   However,  the  density  of  this  plasma  is  very  low 
so  that  the  mean  free  path  for  collisions  is  much  longer 
than  the  dimensions  of  the  obstacle  itself,  the  earth's 
magnetosphere.   Therefore  a  classical  collisional  shock 
wave  cannot  exist.   Satellite  observations  showed  the 
presence  of  a  well-defined  shock  transition  of  the  order 
of  100  Km. 

In  the  theory  of  shock  waves,  it  is  important  to  determine 
the  properties  of  the  medium  after  passage  of  the  shock 
in  terms  of  the  unperturbed  quantities  before  the  shock 
such  as  densities,  pressure  and  temperature.   Assuming  that 
the  shock  profile  is  (quasi-) stationary  in  a  certain 
reference  frame,  a  number  of  relations  between  upstream 
and  downstream  quantities  are  given  by  jump  conditions, 
which  derive  directly  from  the  conservation  of  mass,  momen- 
tum and  energy.   In  usual  gas  dynamics  and  MHD-theory, 
these  jump  conditions  are  sufficient  to  determine  the  basic 
properties  of  the  state  behind  the  shock  although  a  theory 
of  the  shock  structure  and  the  dissipation  process  leading 
to  entropy  production  requires  a  kinetic  description. 

In  case  of  a  collision-free  shock  wave  the  plasma  is, 
in  general,  not  in  local  thermal  equilibrium  behind  the 
shock  transition.   Electron  and  ion  pressures  may  be  quite 
different  from  each  other,  depending  on  the  dissipation 
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mechanism.   Here,  the  jump  conditions  do  not  completely 
describe  the  state  behind  the  shock.   A  dissipation  of  the 
shock  structure  itself  being  needed.   This  would,  e.g. 
provide  information  on  the  ion  or  electron  thermalization 
length. 

Plasma  momentum  transfer  calculated  to  be  collision 
free  was  observed  only  when  a  magnetic  field  was  present 
by  Friedman  and  Dawson  [Ref.  11]  during  plasma  gun  studies. 

Dean  et.  al.  [Ref.  8]  performed  an  experiment  which  was 
similar  to  this  present  one.   A  Nd-doped  glass  laser 
(8J,  30  nsec)  was  used.   The  target  was  0.25  mm  Lucite 
(Cj-H-O-)  fiber  suspended  in  an  ambient  nitrogen  gas.   They 
found  coliisionless  interactions  even  in  the  absence  of 
an  applied  field  and  explained  this  apparent  discrepancy 
by  the  fact  that  spontaneously  generated  fields  always 
existed  in  their  experiment.   They  measured  the  thickness 
of  the  interaction  front  as  approximately  1  mm  and  esti- 
mated the  binary-collision  momentum  transfer  mean  free 
path,  from  a  screened  coulomb  potential  of  the  order  of  a 
centimeter. 

Shock  thicknesses  in  collisional  plasmas  are  usually 
taken  as  one-to-several  collisional  mean  free  paths.   The 
difficulty  of  establishing  a  given  experimental  situation 
as  collision  free  has  been  emphasized  by  Wright  [Ref.  32] 
who  suggested  that  the  results  of  [Ref.  8]  can  be  explained 
by  collisional  coupling  alone.   Dean  et.  al.  [Ref.  9], 
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after  Wright's  suggestion  reperformed  their  experiments. 
They  changed  the  target  to  pure  carbon  filaments  to  find 
out  whether  the  multicomponent  nature  of  the  former  target 
(CgHg02)  was  playing  a  role.   As  a  result  they  found  the  H 
and  0  in  the  target  was  inconsequential.   They  also  changed 
background  gas  and  used  Helium  instead  of  Nitrogen.   The 
basic  phenomena  were  qualitatively  unchanged.   They  sug- 
gested that  the  momentum  coupling  region  which  was  observed, 
was  not  a  shock.   It  was  merely  a  region  of  momentum 
coupling  between  interstreaming  plasmas.   The  mean  free 
path  is  the  e-folding  distance  for  momentum  coupling,  so  that 
several  mean  free  paths  are  required  for  a  high  degree  of  coupling 

The   expanding   plasma   heats  and   transfers  momentum 
to   the   photoionized   background   plasma.    These  inter- 
actions can  be  either  collisional  or  collision-free  depending 
on  densities  and  velocities.   During  an  experiment  per- 
formed by  Stamper  et.  al.  [Ref.  27],  the  background  densi- 
ties have  typically  been  such  that  electron  processes  are 
collisional  while  the  ion  processes  are  collision-free.   They 
observed  a  sharply  defined  momentum  coupling  interaction 
region  between  the  laser  plasma  and  background  plasma,  but 
they  have  not  observed  the  existence  of  a  shock. 

Koopman  and  Tidman  [Ref.  17]  observed  a  thin  collision- 
less  electrostatic  shock  front  propagating  radially  at  the 
leading  edge  of  the  plasma  ball.   In  their  investigation, 
granules  of  LiD  have  been  irradiated  by  15-nsec,  3-J  pulses 
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from  a  Ruby  laser.   The  target  consisted  of  a  cluster  of 
5-to-25  ym  diameter  solid  particles  which  were  injected 
into  the  400-ym  diameter  laser  focal  region  by  an  electro- 
mechanical apparatus.   No  magnetic  fields  were  present. 
The  phenomenon  is  analogous  to  a  blast  wave;  the  expanding 
plasma  acts  like  a  piston  and  sweeps  up  the  background 
plasma.   This  swept-up  plasma  gives  rise  to  a  spherical 
ion-wave  shock. 

The  distinction  between  collision-free  and  collisional 
mechanisms  in  counterstreaming  plasma  studies,  and  the 
appropriate  techniques  for  calculating  collisional  cross 
sections  for  momentum  transfer,  have  been  the  subject  of 
considerable  controversy.   The  question  of  the  role  of 
binary  collisions  is  always  a  very  complicated  one  in  plasma 
experiments.   There  is  still  a  need  for  more  sophisticated 
treatments  of  collisional  effects  in  plasmas. 
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IV.   EXPERIMENTAL  METHODS 

A.  LASER 

The  laser  used  in  this  investigation  was  a  Korad-K-1500 
Q-switched  neodynium-doped  glass  laser.   The  laser  consisted 
of  an  oscillator  and  a  single  amplifier  stage.   The  laser 
was  Q-switched  with  a  Pockels  cell.   The  laser  system  was 
operated  at  a  power  level  of  230  MW,  which  was  obtained 
with  a  pulse  width  (full  width  at  half-maximum)  of  22  nano- 
seconds.  The  output  of  the  system  varied  within  about  10 
percent  of  230  MW  and  no  data  were  accepted  if  obtained  at 
power  levels  outside  this  range.   Other  details  of  the  laser 
are  given  by  Davis  [Ref.  6].   The  laser  focal  spot  size  is 
approximately  0.15  cm  in  diameter  after  focusing  with  a 
30  cm  focal  length  lens. 

B.  LASER  MONITORING  TECHNIQUES 

The  laser  beam  was  monitored  with  a  Korad-K-Dl  photo- 
diode.   This  photodiode  provides  a  signal  which  is  propor- 
tional to  the  laser  power  and  also  by  integrating,  gives 
another  signal  which  is  proportional  to  the  laser  pulse 
energy. 

About  four  percent  of  the  beam  was  reflected  by  a  beam 
splitter  to  a  MgO  diffusing  block  for  viewing  by  the  photo- 
diode.   The  scattered  light  then  passed  through  a  0.1  per- 
cent neutral  density  filter  before  entering  the  photodiode 
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as  shown  in  Figure  2.   The  power  signal  was  used  to  trigger 
a  Tektronix  555  Dual-beam  oscilloscope  and  was  also  dis- 
played on  that  oscilloscope. 

The  energy  signal  was  displayed  on  a  Tektronix  564B 
storage  oscilloscope. 

The  photodiode  energy  signal  was  calibrated  by  L.L. 
McKee  [Ref.  18]  by  using  a  Westinghouse  RN-1  laser  radiometer. 
The  radiometer  provided  an  absolute  measure  of  the  energy 
incident  at  the  target  which  allowed  calibration  of  the 
diode's  integrated  signal. 

The  energy  output  of  the  laser  was  monitored  on  every 
shot.   For  energy  outputs  within  10  percent  of  5.0  Joule, 
the  pulse  width  was  consistently  22  ±  1  nsec.   A  two  minutes 
cooling  period  between  shots  resulted  in  good  reproducibility 

C.   TEMPORAL  REFERENCE  FRAME 

The  temporal  relation  between  the  diagnostic  probe 
response  and  the  leading  edge  of  the  laser  pulse  was  deter- 
mined by  taking  into  account  the  various  cable  delays  and 
laser  beam  optical  path  lengths.   A  long  (10  meter)  delay 
line  was  used  from  the  diagnostic  probe  to  the  oscillo- 
scope in  order  to  separate  the  laser  power  signal  and  early 
diagnostic  responses  which  would  otherwise  be  superimposed 
on  each  other . 

By  assuming  that  the  signal  travels  at  1.98x10    cm/sec 
in  the  RG-174/U  coaxial  cable,  the  delay  was  computed  as 
34  nsec. 
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D.   PLASMA  CHAMBER 

The  target  was  located  in  a  chamber  which  was  specially 
constructed  by  McKee  to  be  used  for  the  investigation  of 
laser-plasma.   The  inside  volume  of  the  chamber  was  cylin- 
drical with  a  diameter  of  10  inches  and  a  height  of  6 
inches  (Figure  3) . 

The  beam  entered  the  chamber  after  passing  through  a 
30  cm  focal  length  lens  and  hit  the  target  at  an  angle  of 
30°  so  that  one  can  probe  the  plasma  along  the  target  nor- 
mal which  appears  to  be  the  axis  of  greatest  symmetry  of 
the  expanding  plasma. 

The  target  used  was  an  aluminum  flat  disc,  1/4  inches 
thick,  and  about  2  inches  in  diameter.   It  was  rotated  after 
a-proximately  10  shots  to  prevent  severe  cratering  of  the 
target  in  the  laser  impact  area. 

The  chamber  was  pumped  to  a  vacuum  of  about  2  x  10    Torr 
of  air  and  then  background  gas  was  allowed  in  by  a  needle 

valve  to  fill  the  chamber  up  to  the  desired  pressure.   The 

—  fi 
vacuum  of  about  1.3  x  10    Torr  was  obtained  by  using  an  oil 

diffusion  pump.   Vacuum  pressures  lower  than  0.1m  Torr 

were  measured  with  an  ionization  gauge.   A  thermocouple 

gauge  was  used  to  measure  pressures  between  0.1m  Torr  and 

1000  m  Torr,  while,  above  1  Torr  pressures  were  monitored 

by  a  manometer  filled  with  diffusion  pump  oil. 
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E.   MAGNETIC  PROBES 

1.   Probe  Calibration 

The  magnetic  probe  used  in  this  investigation  con- 
sisted of  five  turns  of  #40  Formvar  copper  wire.   The  probe 
coil  was  about  1  mm  in  diameter.   The  probe  circuit  response 

was  linear,  with  an  effective  area  (nA)  of  the  inductive 

—  ft      "? 
loop  of  5.175x10    m   for  frequencies  below  about  8  MHz. 

The  effective  area  (nA)  of  the  probe  was  measured  by  using 
Helmholtz  coils.   The  instruments  employed  for  this  measure- 
ment were  a  General  Radio  1330A  continuous  wave  oscillator, 
Hewlett-Packard  411A  r.f.  VTVM,  Helmholtz  coils  and  a 
Tektronix  545B  oscilloscope.   A  circuit  diagram  for  probe 
calibration  is  shown  in  Figure  4. 

The  details  of  construction  of  Helmholtz  coils  and 
the  use  of  these  coils  for  calibration  of  a  probe  is  des- 
cribed by  McLaughlin  [Ref.  19].   The  effective  probe  area, 
according  to  McLaughlin  is  given  by  the  equation 


(HA)   =   k  J-  (1) 

L 


where  K  is  a  parameter  which  differs  only  from  coil  to 

coil  (for  the  coil  which  was  used  in  this  investigation, 

2 
K  =  0.025  m  ) ,  E  is  the  voltage  across  the  probe  connection 

cable  and  V   is  the  voltage  applied  to  the  Helmholtz  coil. 

2.   Probe  Signal  Reliability  and  Sensitivity 

To  insure  that  the  magnetic  probe  signals  (9B/3t) 

were  meaningful,  several  steps  were  taken.   The  probe 
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signals  were  monitored  under  experimental  conditions  for 
which  no  target  was  in  place,  to  verify  that  a  null  signal 
was  obtained.   Also,  the  probe  was  rotated  by  180°  about 
its  axis  and  the  probe  signals  were  checked  to  make  sure 
that  they  reverted  in  polarity.   This  demonstrated  that 
the  signals  were  not  electrostatic  in  nature.   Magnetic 
field  signals  were  recorded  on  a  Tektronix  555  dual-beam 
oscilloscope. 

3 .   Probe  Signal  Integration 

According  to  electromagnetic  theory,  the  output  of 
a  coil  is 


e(t)   =   nA  ||  (26 


where  r|A  is  the  effective  coil  area,  B  is  the  magnetic 
field,  and  e(t)  is  the  voltage  generated  across  the  coil 
which  was  recorded  on  a  Tektronix  555  dual-beam  oscillo- 
scope.  Integration  of  equation  (26)  gives: 


B  t 

B   =    /  dB   -   i   /  e(t)dt  (27) 

0  0 


t 
The  integration  of  the  probe  signal  (  /  e(t)dt)  was 

0 
performed  by  measuring  the  area  under  the  e(t)  curve. 
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F.   ELECTRIC  DOUBLE  PROBES 

Double  probe  thoery  can  be  found  in  several  references 
[Ref.  3]  and  [Ref.  4]. 

At  sufficiently  high  bias  voltage,  the  probe  draws 
the  ion  saturation  current  I 


An  n,  Ze   V_  (28 

1   1       Lp 


where  A,  is  the  probe  area  perpendicular  to  the  direction 


of  plasma  flow,  n  •  Ze  is  the  average  ionic  charge  per  unit 
volume  in  the  laser  plasma  and  V   is  the  magnitude  of  the 
plasma  flow  velocity. 

In  equation  (28)  the  assumptions  are: 

(1)  The  areas  of  probe  tips  are  equal. 

« 

(2)  The  plasma  parameters  are  uniform  across  the  probe 
tips. 

(3)  Plasma-probe  interactions  such  as  secondary  elec- 
tron emission  can  be  neglected. 

The  secondary  electron  emission  coefficient  (e)  due  to 
ion  bombardment  (y.),  defined  as  the  average  number  of 
electrons  ejected  per  incident  ion,  has  been  measured  by 
Koopman  [Ref.  16]  for  similar  experimental  conditions  and 
found  to  fall  in  the  range  of  0.0  -  0.5.   Therefore  secon- 
dary electron  emission  will  contribute  to  the  total  current 
drawn  between  the  probe  tips.   After  taking  secondary 
electron  emission  into  account  equation  (28)  becomes 
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I  ~  (1+e)  A1  niZe  VLp  (29) 


The  probe  was  biased  at  20  volts  to  collect  a  saturated 
ion  current.   A  large  secondary  emission  coefficient  of 
e  =  0.5  was  used  in  obtaining  a  quantitative  idea  of  plasma 
density. 

The  probe  circuit  is  shown  in  Figure  (5).   The  probe 
circuitry  was  enclosed  by  an  aluminum  housing  attached  to 
the  probe  and  the  entire  system  was  isolated  from  ground. 
The  diameter  of  the  cylindrical  tips  was  about  1  mil,  they 
had  a  length  of  1  mm  and  were  about  1  mm  apart. 

The  double  probe  signals  were  recorded  on  a  Tektronix 
555  dual-beam  oscilloscope.   To  obtain  an  estimate  of  the 
plasma  density,  equation  (t)  is  solved  for  n •  and  evaluated 
with  Z  =  1.5,  £  =  0.5,  A,  =  2.54 x  10~   cm2 ,  I  =  ^  ,  R  =  1  Q 


,  nn    in22      V  (volts)  ,-,,., 

n  .   =   1.09  x  10    x  - t t —, —        30 

'i  V     cm/sec 


where  V  is  the  voltage  across  a  1  ohm-resistor  (see  Figure  5) 

G.   DATA  COLLECTION  AND  ERROR  ESTIMATES 

The  magnetic  probe  was  inserted  into  the  vacuum  chamber 
through  a  small  pipe  which  was  soldered  onto  a  brass  base 
plate. 

The  base  plate  could  be  moved  about  on  an  "0"  ring  seal 
located  on  the  top  of  the  chamber  without  breaking  the 
vacuum  seal.   The  base  plate  was  scribed  with  reference  lines 
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These  lines  allowed  positioning  of  the  plate  over  a  one 
millimeter  grid  attached  to  the  top  of  the  chamber.   The 
position  of  the  probe  relative  to  the  target  was  measured 
internally  at  the  beginning  of  each  data  session  to  insure 
that  the  probe  had  not  moved  in  its  holder.   This  method 
allowed  positioning  of  the  probe  to  within  ±  -j  mm  of  the 
stated  position  in  relation  to  the  target.   The  electric 
probe  was  inserted  into  the  chamber  from  a  side  window. 

The  relative  measure  of  plasma  density  from  point  to 
point  is  considered  good  due  to  the  high  reproducibility 
of  the  laser  pulses.   It  is  estimated  that  these  relative 
values  are  accurate  to  within  ±  10  percent. 

The  laser  gave  very  reproducible  pulses;  the  output 
energy  was  allowed  to  vary  as  much  as  ±  10  percent. 

The  area  of  probe  tips  is  considered  accurate  to  within 
10  percent.   The  plasma  expansion  velocity  is  considered 
accurate  to  within  10  percent. 
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V.   EXPERIMENTAL  RESULTS 

During  this  investigation,  hydrogen  gas  was  used  as 
the  background  gas.   An  electric  probe  was  introduced  from 
the  port  No.  6.  (Figure  3)  and  the  magnetic  probe  from  the 
top  window. 

To  decelerate  the  plasma  front  rapidly ,  a  glass  plate 
is  placed  into  the  vacuum  chamber  (Figure  7) .   The  main 
influence  of  the  glass  plate  on  spontaneously  produced 
magnetic  field  is  displayed  in  Figure  28.   This  figure 
indicates  that  the  maximum  value  of  reverse  field  when  the 
glass  plate  is  in  place  (at  z  =  1.15  cm) ,  measured  with 
the  probe  at  r  =  0.4  cm,  0=0°  and  z  =  1.0  cm,  is  greater 
than  the  maximum  value  of  reverse  field  when  the  glass 
plate  is  not  in  place  (at  z  =  7.5  cm) . 

Figures  8  and  9  display  the  relation  of  3B/3t  to  the 

—  ft 
plasma  density  for  a  vacuum  pressure  of  2  x  10    Torr,  as 

a  function  of  time.   When  the  glass  plate  is  in  place,  the 

magnitude  of  3B/3t  is  greater  than  the  magnitude  of  3B/3t 

when  the  glass  plate  is  not  in  place. 

Figure  16  demonstrates  that  the  slope  of  the  plasma 
density  curve  when  the  glass  plate  is  at  z  =  1.15  cm  is  not 
very  different  from  the  slope  of  the  plasma  density  curve 
when  the  glass  plate  is  at  z  =  7 . 5  cm  . 

Figure  17  shows  that  when  the  background  gas  pressure 
is  0.5  Torr  (H2) ,  the  existence  of  the  glass  plate  at 
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z  =  1.15  cm  causes  a  noticeable  steepening  and  a  sudden 
increase  in  plasma  density. 

Figures  10  and  11  show  the  relation  of  9B/8t  to  the 
plasma  density  at  a  background  pressure  of  0.5  Torr  (H^) 
as  a  function  of  time  for  two  locations  of  the  glass  plate. 

When  the  background  gas  pressure  is  5.0  Torr  (H2),  the 
sudden  increase  in  plasma  density  when  the  glass  plate  is 
in  place  (Figure  18)  is  much  more  pronounced  (see  also 
Figures  13,  15  and  28).   For  5.0  Torr  (H2)  background 
pressure,  the  difference  between  the  maximum  value  of 
reverse  field  when  the  glass  plate  is  at  z  =  1.15  cm  and  the 
maximum  value  of  reverse  field  when  the  glass  plate  is 
at  z  =  7.5  cm  is  greater  than  the  difference  between  the 
same  quantities  for  a  background  pressure  of  0.5  Torr  (EU) . 

Figures  12  and  14  display  the  relation  of  3B/3t  to  the 
plasma  density  for  5.0  Torr  background  gas  pressure  at  the 
point  r  =  0.4  cm,  6=0°  and  z  =  1.0  cm  for  two  different 
locational  of  the  glass  plate  (note  that  each  curve  in  this 
figure  was  multiplied  by  a  different  factor) . 

Figures  19,  21  and  23  demonstrate  the  comparison  of  two 
plasma  density  curves  each  for  different  background  gas 
pressure,  at  the  point  r  =  0.4  cm,  9=0°  and  z  =  1.0  cm 
with  the  glass  plate  at  z  =  1.15  cm.   The  slope  of  each  of 
these  two  curves  is  almost  the  same.   The  comparisons  of 
the  same  curves  when  the  glass  plate  is  at  z  =  7.5  cm  are 
displayed  in  Figures  20,  22  and  24.   From  these  figures,  one 
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can  notice  that  the  slopes  of  the  two  density  curves  each 
one  for  different  background  gas  pressures  are  different 
from  each  other. 

Figures  25,  26  and  27  illustrate  the  comparison  of 
the  magnetic  probe  signals  (3B/3t)  when  the  glass  plate 
is  located  at  z  =  1.15  cm  with  the  magnetic  probe  signals 

(3B/3t)  when  the  glass  plate  is  located  at  z  =  7.5  cm  for 

—  ft 
2x10    Torr,  0.5  Torr  (H2)  and  5.0  Torr  (H2)  background 

pressures.   From  these  figures,  one  can  see  that  the  reverse 
field  maximum  is  the  greatest  when  the  background  gas 
pressure  is  5.0  Torr  (H2)  and  glass  plate  is  at  z  =  1.15  cm. 

Figures  29,  31  and  33  illustrate  the  comparison  of  the 
two  magnetic  probe  signals  (3B/3t),  each  one  for  a  differ- 
ent background  pressure,  as  a  function  of  time  at  the  point 

« 
r  =  0.4  cm,  8=0°  and  z  =  1.0  cm  and  the  glass  plate  is 

in  place.   An  investigation  of  these  figures  indicates  that 
the  higher  the  background  gas  pressure,  the  greater  the 
area  under  the  negative  portion  of  the  3B/3t  curve,  while 
the  area  under  the  positive  portion  of  3B/3t  curve  becomes 
smaller.   This  means  that  as  the  background  gas  pressure 
is  increased,  the  reverse  field  maximum  becomes  greater  when 
the  glass  plate  is  in  place.   Figures  30,  32  and  34  illus- 
trate the  same  comparison  when  the  glass  plate  is  at 
z  =  7.5  cm. 

For  a  vacuum  pressure  of  2x10    Torr,  the  magnitude 
of  the  magnetic  probe  signal  (3B/3t)  when  the  glass  plate 
is  in  place  is  about  4  times  greater  than  the  magnitude  of 
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the  magnetic  probe  signal  (3B/3t)  when  the  glass  plate  is 
not  in  place  (see  Figure  25) .   (Note  that  each  curve  in 
this  figure  was  multiplied  by  a  different  factor.) 

When  the  glass  plate  is  in  place,  the  magnitude  of  the 
magnetic  probe  signal  (3B/3t)  at  a  vacuum  pressure  of  2x10 
Torr  is  greater  than  the  magnitude  of  the  magnetic  probe 
signals  (3B/3t)  at  0.5  Torr  (H2)  background  pressure  and 
at  5.0  Torr  (H2)  background  pressure  (see  Figures  29  and 
31).   But  when  the  glass  plate  is  not  in  place,  the  magnitude 
of  the  magnetic  probe  signal  (3B/3t)  at  2x10    Torr  vacuum 
pressure  is  smaller  than  the  magnitude  of  the  magnetic 
probe  signals  (3B/3t)  at  0.5  Torr  (H2)  background  pressure 
and  at  5.0  Torr  (H2)  background  pressure  (see  Figures  30 
and  32).   From  Figure  34,  it  can  be  seen  that  when  the  glass 
plate  is  not  in  place,  the  magnitude  of  the  magnetic  probe 
signal  (3B/3t)  at  0.5  Torr  (JU)  background  pressure  is 
greater  than  the  magnitude  of  the  magnetic  probe  signal 
(3B/3t)  at  5.0  Torr  (H2)  background  pressure.   (Note  that 
each  curve  in  this  figure  was  multiplied  by  a  different 
factor.)   But  Figure  33  illustrates  that  when  the  glass 
plate  is  in  place  the  magnitude  of  the  magnetic  probe  signal 
(3B/3t)  at  a  background  pressure  of  0.5  Torr  (H2) ,  between 
0  and  100  nsec  is  greater  but  between  100  and  4  50  nsec  it 
is  smaller  than  the  magnitude  of  the  magnetic  probe  signal 
(3B/3t)  at  a  background  pressure  of  5.0  Torr  (H2).   (Note 
that  each  curve  in  this  figure  was  multiplied  by  a  different 
factor . ) 
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Figure  35  illustrates  that  when  the  glass  plate  is 
in  place,  for  a  background  pressure  of  5.0  Torr  (H~),  the 
slope  of  the  plasma  density  curve  at  the  point  r  =  0.4  cm, 
9=0°  and  z  =  1.0  cm  is  greater  than  the  slope  of  the 
plasma  density  curve  at  the  point  r  =  0.8  cm,  6=0°  and 
z  =  1.0  cm.   But  when  the  glass  plate  is  not  in  place 
Figure  36)  the  slopes  of  two  plasma  density  curves  are  the 
same . 

During  this  investigation,  for  a  background  pressure  of 
5.0  Torr  (EU) /  the  maximum  value  of  the  magnitude  of  the 
reverse  magnetic  field  was  measured  as  150  gauss  when  the 
glass  plate  is  in  place  at  the  point  r  =  0.4  cm,  9=0° 
and  z  =  1.0  cm  and  at  t  =  440  nsec ,  while  at  same  background 
pressure,  point  and  time  but  when  the  glass  plate  is  not 
in  place,  the  maximum  value  of  the  magnitude  of  the  reverse 
magnetic  field  was  measured  as  45  gauss  only. 
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VI.   COMPARISON  OF  THEORY  AND  EXPERIMENTAL  RESULTS 

As  mentioned  in  Section  III-C,  the  production  of  initial 
field  requires  a  density  gradient  in  the  -z  direction  and 
a  temperature  gradient  in  the  -r  direction  (see  Figure  6)  . 
The  cross  product  of  these  two  vectors  generates  a  magnetic 
field  in  the  6  direction  i.e.,  in  the  clockwise  direction 
when  looking  in  the  direction  of  the  laser  beam. 

If  one  of  the  gradients  is  reversed  to  the  opposite 
direction,  a  magnetic  field,  in  the  counterclockwise 
direction  is  produced.   The  radial  temperature  gradient  is 
always  in  the  -r  direction  and  its  direction  remains  unchanged 
either  the  glass  plate  is  in  place  or  not.   Therefore  the 
only  possible  gradient  which  can  be  reversed  by  placing 
the  glass  plate  at  z  =  1.15  cm,  is  the  axial  density  gradi- 
ent.  The  equation  of  continuity  is 


|I1  +  V  •  (nv)   =  o  (31) 


if  cylindrical  symmetry  is  assumed; 


8n    3n  .  dz  .  3r»  .  dr 

3t    3z    dt    3r    dt 


0  (32 


Also  for  the  laser-produced  plasma  expanding  in  the  z- 
direction 


¥■  »   P- 

dt       dt 


thus:   -r—  ~  ~  TT  *  v  where  V  is  the  expansion  velocity  of 

the  plasma  front. 

Figure  13  demonstrates  the  relation  of  Bq  to  the 

plasma  density  as  a  function  of  time  at  the  point  r  =  0.4  cm, 

0=0°  and  z  =  1.0  cm  for  a  background  gas  pressure  of 

5.0  Torr  (H2)  when  the  glass  plate  is  at  z  =  1.15  cm. 

Between  t  a  100  nsec  and  t  z    180  nsec.  ~   is  positive, 

d  t 

which  means  that  |p-  is  negative  or  the  axial  density 
gradient  is  in  the  -z  direction.   Figure  28  illustrates 
that  at  t  ?  150  nsec  the  reverse  field  appears  when  the 
glass  plate  is  in  place.   Since  the  axial  density  gradient 
at  this  time  is  in  the -z  direction,  the  production  of  the 
reverse  field  can  not  be  correlated  to  any  change  of  sign 
in  the  direction  of  the  axial  density  gradient,  since 
there  is  no  such  change  when  the  reverse  field  is  "turned 
on . 

There  always  exists  a  radial  density  gradient  in  the  -r 
direction,  no  matter  where  the  glass  plate  is  (see  Figures 
35  and  36).   Therefore  it  requires  an  axial  temperature 
gradient  in  the  -z  direction.   Such  an  axial  temperature 
gradient  can  be  generated  by  the  shock  heating  of  the 
plasma. 

Figure  2  8  demonstrates  that  the  magnitude  of  the 
initial  magnetic  field,  when  the  glass  plate  is  in  place 
is  smaller  than  the  magnitude  of  the  initial  magnetic 
field  when  the  glass  plate  is  not  in  place.   This  decrease 
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in  magnitude  of  the  initial  magnetic  field  is  caused  by 
the  decrease  in  axial  density  gradient  due  to  pile-up  of 
laser  plasma  in  front  of  the  glass  plate.   This  piled 
up  plasma  then  diffuses  radially  and  also  causes  a  decrease 
in  radial  density  gradient.   This  decrease  in  radial  den- 
sity gradient  can  be  seen  if  one  compares  Figure  35  and 
Figure  36. 

Figure  18  displays  the  steepening  in  the  plasma  den- 
sity between  150  and  2  00  nsec,  when  the  expanding  plasma 
impinges  on  a  glass  plate.   This  steepening  can  be  related 
to  the  formation  of  a  shock  wave.   Such  a  shock  wave  causes 
an  increase  in  plasma  temperature  (see  Appendix  A)  and  this 
increase  results  in  an  axial  temperature  gradient  in  the 
-2  direction,  since  the  shock  wave  travels  in  the  +z  direction 

As  long  as  the  ratio, 

i  *Bi   i 

7  T  I  I  V  n  I 


9t      '  r  e M  z 


dB  V  T    V  n 

I r  I    '  z  e  M  r  ' ' 

1  at  ' 

is  greater  than  one,  the  reverse  field  only  reduces  the 
magnitude  of  the  initial  magnetic  field.   But  when  this 
ratio  decreases  to  a  value  which  is  less  than  one,  the 
reverse  field  becomes  the  dominant  field. 

From  Figure  18,  the  ratio  of  two  densities,  one  at 
t  =  159.37  nsec  the  other  at  t  =  180.83  nsec,  is: 


—  =   1.5      and      At   =   21.46  nsec 

ni 
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in  order  to  obtain  this  density  ratio,  the  shock  wave 
must  have  a  Mach  number  of  the  order  of  1.37.   According 
to  the  classification  of  Jaffrin  and  Probstein  [Ref.  15] 
this  shock  is  a  "moderate"  shock  whose  thickness  is  of 

the  order  of  H/e,    where  I    is  the  downstream  (before  shock) 

me  1/2 
ion-ion  collision  mean  free  path,  and  e  =  ( — )    .   The 

i 

electron  temperature  increases  in  this  layer  but  the  ion 

temperature  increases  in  the  second  inner  layer  (see  Section 
III-E)  whose  approximate  thickness  is  the  downstream  ion- 
ion  collision  mean  free  path  I . 

For  this  Mach  number  M  =  1.37,  a  temperature  ratio, 


el     ll 


is  obtained  from  the  Rankine-Hugoniot  jump  conditions 
(see  Appendix  A,  equation  3) . 

Since  the  shock  thickness  can  be  defined  from  the 
electric  probe  data  as  5  =  Vt ,  where  V  is  the  velocity  of 
the  expanding  front  and  t  is  the  rise  time  of  the  probe 
signal. 

If  no  interactions  were  occurring  between  the  inter- 
streaming  plasmas,  then  the  plasma  expansion  velocity 
would  be  independent  of  background  gas  pressure.   The 
observed  dependence  of  expansion  velocity  can  be  seen  in 
Figures  20,  22,  24,  between  100  and  200  nsec.   From  these 
figures  one  can  conclude  that  as  background  gas  pressure 
increases,  the  expansion  velocity  decreases. 
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To  have  an  idea  about  the  structure  of  the  shock 
one  must  compare  the  shock  thickness  with  the  collison 
mean  free  path  length.   The  ion-ion  collision  mean  free 
path  can  be  calculated  from 


1013  2 

I    (cm)  z      =^-= x  T.  (eV)  , 

n  •  (cm   )  x  in  a 
1 

where  T.  is  the  ion  temperature.   In  order  to  calculate  I, 

one  needs  to  know  n •  and  T . . 

'i       1 

1    2 
From  kT.  =  -^m . V  ,  where  V  is  the  expansion  velocity 

and  m.  is  the  ion  mass,  T.  can  be  estimated.   If  one 

assumes  that  the  Shock  wave  propagates  in  the  background 

gas,  or  background  gas  flows  into  the  shock  relative  to 

the  frame  of  the  shock,  for  hydrogen  background  gas, 

-24  7 

m.  =  m   =  1.673  x  10     g  and  V  z    10   cm/sec.   For  these 

values  of  m.  and  V,  T.  z    52  eV  is  obtained.   Assuming  a 

fully  dissociated  and  fully  ionized  background  gas  at 

5.0  Torr  the  downstream  ion  density  is 


n-      2x3.3x5x10    cm        3.3x10    cm 


For  these  values  of  n •  and  T. ,  i    is  of  the  order  of 

i      l 

8.2  x  io"3  cm. 

The  shock  thickness  can  be  obtained  from  Figure  18, 

-9  7 

in  which  t  =  At  =  21.46  x  10    sec  and  for  V  ~  10   cm/sec 

-2 
6  =  xV  =  21.46  x  io    =0.21  cm.   Since  the  shock  thickness  6 
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is  greater  than  the  ion-ion  collision  mean  free  path, 
the  formation  of  a  collisional  shock  wave  must  be  expected 
In  summary  it  can  be  said  that  whenever  the  background 
gas  pressure  is  high  enough,  the  expanding  laser  plasma 
acts  like  a  piston  and  sweeps  up  the  photoionized  back- 
ground plasma.   In  front  of  this  piston  a  shock  is  formed 
to  drive  an  axial  temperature  gradient  which  is  large 
enough  to  produce  a  reverse  magnetic  field. 
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VII.   SUMMARY  AND  CONCLUSIONS 

The  observed  dependence  of  the  plasma  front  expansion 
velocity  on  background  gas  pressure  indicates  the  exis- 
tence of  an  interaction  between  the  laser  plasma  and  the 
photoionized  background  plasma.   As  a  result  of  this  inter- 
action, a  collisional  shock  wave,  driven  by  the  laser 
plasma  piston  is  formed.   This  shock  wave  heats  and  com- 
presses the  plasma. 

Initially  the  spontaneously  produced  magnetic  field  is 
in  the  6-direction  (see  Figure  1)  and  the  thermal  source 
term  for  these  magnetic  fields  consists  of  the  cross 
product  of  the  radial  density  gradient  and  the  axial  den- 
sity gradient.   When  the  shock  wave  is  formed  an  axial 
temperature  gradient  is  produced.   This  axial  temperature 
gradient  and  the  radial  density  gradient  which  always 
exists  in  the  plume  shaped  laser  plasma  are  the  main  cause 
of  the  reversal  of  these  initial  magnetic  fields. 

The  comparison  of  the  shock  thickness  and  the  colli- 
sion mean  free  path  determines  the  type  of  the  shock.   If 
the  shock  thickness  is  greater  than  the  collision  mean 
free  path,  the  shock  is  collisional  in  nature. 

Determining  the  mean  free  path  length  in  a  plasma  is 
always  a  complicated  undertaking  and  the  results  should  be 
considered  as  crude  estimates.   More  sophisticated  treat- 
ments are  necessary  in  order  to  determine  the  collisional 
effects  in  plasmas. 
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The  reversal  of  the  spontaneous  magnetic  field  indi- 
cates the  existence  of  a  shock  wave  propagating  through 
the  plasma.   Such  generation  of  magnetic  fields  by  shock 
waves  can  be  a  very  important  mechanism  in  astrophysics. 
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VIII.   RECOMMENDATION  FOR  FURTHER  STUDY 

The  electron  temperature  has  been  evaluated  indirectly 
from  the  shock  parameters.   An  experimental  measurement 
of  the  electron  temperature  gradient  is  required  for  a 
further  study  of  field  reversal  in  shock  waves.   Plasma 
temperature,  before  and  after  the  plasma  front  encounters 
the  glass  plate  can  be  measured  by  means  of  x-ray  diagnostic 
techniques  or  by  other  means.   Also,  more  accurate  data 
about  the  r  dependence  of  the  plasma  density  are  needed 
to  provide  a  quantitative  comparison  between  experiment 
and  theory.   The  data  about  the  r  dependence  of  the  plasma 
density  can  be  obtained  by  use  of  the  electric  probe 
technique. 
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APPENDIX  A 

Rankine-Hugoniot  jump  conditions  can  be  determined  by 

using  the  equation  of  continuity,  conservation  of  momentum 

density  and  conservation  of  energy  density  in  the  plasma 

for  steady  state  conditions  [Ref.  29]. 

c        5 

The  Mach  number  is  defined  for  — *-  =  y  =  -=-  as : 

CV        3 

vx 

M   =    5k(T  1+T.n)  (1) 

r el    11   -j 

L     3M      J 

where  V,  is  the  downstream  (before  shock)  fluid  velcoity, 
T  .  is  the  downstream  electron  temperature,  T? ,  is  the 
downstream  ion  temperature  and  M  is  the  ion  mass. 
The  jump  conditions  are: 


nl      V2      13 
n2      vl      4  Mz 


T   +T  2        2 

i2   e2  _   (MN-3)  (5M  -1)  ... 

t.-+t  :  "   T^p. (3) 

ll   el  16M 


where  the  subscript  "2"  implies  upstream  (after  shock) 
value. 
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FIGURES 
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Figure  1.   Cylindrical  polar  coordinate  system 
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Figure  2.   Block  diagram  of  experimental  layout 
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Figure  3.   Top  view  of  vacuum  chamber.   Port  #1  is  the  laser 
beam  entry  port,  #2  is  the  reflected  laser  beam 
observation  port,  #3  is  the  transmitted  laser 
beam  observation  port,  #4  is  the  optical  observa- 
tion port,  #5  is  the  glass  plate  entry  port, 
#6  is  the  electric  probe  entry  port,  #7  is  the 
probe  observation  port. 
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Figure  5 .   Electric  double  probe  circuit 
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Figure  6 .   Geometry  of  a  laser  produced  plasma  and 
the  self  generated  magnetic  field 
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FIGURE  7.   Arrangement  of  the  target,  glass  plate 
and  probes 
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